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A  New  Method  of  Measuring  Multigrade  Oil 

Shear  Elasticity  and  Viscosity  at 

High  Shear  Rates 


ABSTRACT 

Recent  journal  bearing  oil  film  thickness 
data  (1)*  strongly  suggest  that,  in  order  to 
develop  ■ininuB-viscoaity  aultigrade  oils  for 
■axinizing  fuel  econoay,  it  is  not  enough  to 
exaaine  viscosity  n  alone,  even  at  150C  at  a 
shear  rate  f  =  10*a~>:  it  is  essential  also  to 
exaaine  a  certain  kind  of  shear  elasticity,  of 
which  one  aeaaure  is  the  first  normal  atresa 
difference  Ni  in  shear  flow.  Within  a  group  of 
isoviscoua  aultigrade  oils,  Ni  values  can  differ 
significantly;  for  lOW-40  A  15W-40  oils,  an 
average  of  75S  of  the  ainiaua  oil  fila  thickness 
can  be  attributed  to  Ni  (1).  Here  we  describe 
(a)  the  aethod  used  in  (1)  to  obtain  high-f  Ni 
data;  (b)  an  increase  in  the  Ni -aeasureaent 
upper  liait  faax  froa  3  10®s'i  to  10®s-‘  (where 
n  =  lOcP);  &  (c)  an  increase  in  the  n~weasure- 
aent  Tsax  to  5  10*s~^.  This  is  close  to  the  Taax 
value  found  in  (1)  A  is  5  tiaes  greater  than 
that  at  present  attainable  with  other  cawKr- 
cially  available  viscoaeters. 


IT  IS  WBU  nOWN  THAT,  in  the  first  instance, 
polyaers  were  added  to  lubricating  oils  in  order 
to  reduce  the  teaperature  dependence  of 
viscosity.  However,  the  addition  of  polyaer  has 
other  effects,  notably  the  creation  of  a  group 
of  non-Newtonian  rheological  properties  that  can 
collectively  be  termed  elastic,  which  are  not 
exhibited  by  single  grade  oils  (containing  no 
polyaeric  additives).  The  possible  beneficial 
effects  of  oil  elasticity  on  Journal  bearing 
lubrication  have  been  considered  for  over  30 
years  (2-6);  until  recently,  the  experiaental 
investigations  (aostly  of  wear)  have  been 
somewhat  inconclusive. 

'Numbers  in  parentheses  (not  preceded  by  Bq.) 
refer  to  the  RBFBRBNCBS  sectimi  at  the  end  of 
the  paper. 
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FIGURB  1  Addition  of  polyaer  to  lubricating  oil 
generates  an  elastic  contribution  Ni  to  the 
state  of  stress  in  shear  flow  in  full  fila 
hydrodynaaic  lubrication  in  a  journal  bearing. 
The  elastic  contribution  has  the  nature  of  a 
tensile  hoop  stress  acting  tangentially*  to  the 
shear  lines.  Where  the  gap  is  saaller,  the  shear 
rate  A  Ni  values  are  larger.  The  Ni  distribution 
thus  gives  a  resultant  force  on  the  journal, 
which  acts  in  line  with  the  ainiaua  fila 
thickness  disaster  in  a  direction  tending  to 
increase  the  ainiaua  fila  thickness. 

For  aultigrade  oils  in  shear  flow,  the 
state  of  stress  in  a  liquid  element  can  be  re¬ 
garded  (7)  as  a  superposition  of  an  isotropic 
pressure  p,  a  shear  stress  o,  A  first  A  second 
normal  stress  differences  Ni  A  Nz .  The  rate  at 
which  surface  tractions  do  work  on  the  element 
is  equal  to  ot  per  unit  volume  A  does  not  de¬ 
pend  on  Ni .  For  most  non-Newtonian  polymeric 
liquids,  Ni  >  0.  Ni  is  related  to  elastic  re¬ 
covery  at  constant  volume  (8).  Ni*/a,  one  def¬ 
inition  of  an  oil  shear  modulus,  is  of  the  order 
lO^-lO'Pa  A  is  a  quantity  distinct  from  the  bulk 
Bodulus  (typically  of  the  order  10*-10iopa). 

Fig.l  shows  one  possible  aechanisa  whereby, 
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for  the  case  of  full  fil*  hydrodynaaic  lubri¬ 
cation,  oil  elasticity  can,  at  least  in  prin¬ 
ciple,  generate  an  additional  force  on  the 
Journal  in  a  direction  tending  to  increase  the 
■iniaua  oil  film  thickness.  This  elastic  force 
contribution  is  at  right  angles  to  that  due  to 
viscosity  &  is  of  particular  value  because, 
being  associated  solely  with  Ni ,  the  first  nor- 
■al  stress  difference  (a  conservative  contri¬ 
bution  to  stress),  no  energy  is  consuaed  in  its 
generation.  If,  therefore,  significant  elastic 
effects  can  be  generated  fc  aaintained  under 
practical  operating  conditions,  it  should  be 
possible  to  reduce  oil  viscosity  while  aain- 
taining  acceptable  oil  fila  thicknesses.  If 
Journal  bearing  failure  is  a  critical  factor 
where  oil  viscosity  reduction  is  to  be  con¬ 
sidered  in  relation  to  coaplete  engine  perform¬ 
ance,  it  is  reasonable  to  hope  that  acceptable 
lower-viscosity  aultigrade  oils  could  be  formu¬ 
lated  &  hence  that  energy  savings  could  result. 

In  shear  flow  between  vertical  cylinders  in 
relative  rotation,  Ni  acts  like  a  hoop  stress 
giving  rise  to  the  "Weissenberg  Effect"  in  which 
the  liquid  cliabs  up  the  inner  cylinder.  If  the 
cylinders  are  coaxial,  this  hoop  stress  distri¬ 
bution  has  rotational  syHetry  about  the  coaaon 
axis  &  so  gives  no  resultant  force  contribution 
in  planes  normal  to  the  cylinders.  In  a  Journal 
bearing  under  load,  the  Journal  &  bearing  are 
not  coaxial  (fig.l)  &  there  is  a  resultant 
"elastic"  force  contribution  froa  the  asyMetric 
Ni  distribution. 


FIGURE  2  The  above  recent  GM/Shell  data  (I) 
show  that,  in  order  to  correlate  aeasured  val¬ 
ues  of  ainiaua  oil  fila  thickness  in  a  front 
aain  bearing  of  a  V-6  engine  running  at  3000  rpa 
with  rheological  properties  of  the  oils,  it  is 
necessary  to  include  the  elastic  quantity  Ni  as 
well  as  the  viscosity  q.  Oil  viscosities  in  aPas 
(cP)  are  given  on  the  bottom  line.  On  average, 
at  equal  high-teaperature,  high-shear-rate  vis¬ 
cosities,  the  ainiaua  fila  thicknesses  for  the 
SAE  lOlf-40  &  15W-40  elastic  oils  are  75%  greater 
than  those  for  the  SAE  40  inelastic  oil. 

Recently  (1),  electrical  resistance  probes 
were  used  to  aeasure  oil  fila  thickness  in  a  V-6 


front  aain  bearing  at  3000  rpa.  A  concentric  cy¬ 
linder  apparatus  was  used  to  aeasure  viscosity. 
A  new  instniaent,  the  LODGE  STRESSMBTER^ "  >  FOR 
HIGH  SHEAR  RATES,  was  used  to  aeasure  Ni  for 
shear  rates  up  to  3  lO^s'^  &  teapteratures  up  to 
60C;  extrapolation  was  used  in  order  to  extend 
the  range  to  ISOC  It  6  lO^s'^,  the  highest  est¬ 
imated  shear  rate  used.  For  the  set  of  2  single 
grade  &  13  aultigrade  oils  investigated,  in¬ 
clusion  of  an  Ni  tera  in  a  regression  equation 
for  the  ainiaua  oil  fila  thickness  gave  an 
increase  in  the  correlation  coefficient  R^  from 
0.15  to  0.70  compared  with  the  case  in  which 
viscosity  alone  was  used.  Some  results  (froa 
Table  12  of  (1))  are  presented  in  fig. 2.  The 
authors  conclude  that 

"multigrade  oils  reault  in  thicker  oil  filme 
than  mingle  grades  at  equal  high-temperature, 
higb-ahear-rate  viacoai ties; .film  thicknesses 
for  the  SAB  lOh-40  A  lStl-40  oils  are,  on 
average,  7St  greater  than  the  SAB  40  single 
grade  oil."  (our  italics) 

The  importance  of  elastic  contributions 
shown  by  these  results  is  surprising  in  view  of 
theoretical  estimates  (4-6)  which  suggest  that 
values  of  the  stress  ratio  Ni/o  (where  o  denotes 
the  shear  stress)  would  need  to  be  of  the  order 
of  several  hundred  to  sake  the  fig.l  Hi  contri¬ 
bution  comparable  in  magnitude  with  the  viscous 
contribution.  Since  the  measured  values  of  Ni/o 
were  below  10,  the  observed  strong  correlation 
between  fila  thickness  &  Hi  data  is  unexpected 
and  has  yet  to  be  explained. 

In  seeking  an  explanation,  possible  factors 
to  be  considered  include  the  following:  the  aax- 
iaua  value  0.96  of  the  eccentricity  ratio  (9)  in 
the  aeasurements  was  such  greater  than  the  val¬ 
ues  used  in  the  calculations;  squeeze  fila  be¬ 
havior  associated  particularly  with  dynamic 
loading  of  short  bearings  aay  play  a  significant 
role  (10);  polymer  additives  aight  affect  cav¬ 
itation  &  boundary  lubrication;  Hi  nay  play  the 
sole  role  of  a  sensitive  detector  of  high  MH 
polymer  or  of  polyner-polyner  interactions,  & 
the  polymer  itself  aay  act  (in  some  as  yet  un¬ 
known  manner)  to  increase  the  ainiaua  film 
thickness.  Examples  are  known  (11,  12)  in  which 
polymer  MH  changes  affect  Hi  but  not  n* 

The  statistical  significance  of  the  cor¬ 
relation  between  Ni  !■  ainiaua  fila  thickness  is 
independent  of  the  fact  that  we  have  as  yet  no 
adequate  hypothesis  for  a  mechanism  which  aight 
account  for  this  correlation. 

Whatever  the  mechanism  aay  be,  it  seems  to 
be  beyond  reasonable  doubt  that  the  GM/Shell  re¬ 
sults  (1)  establish  two  important  facts  in  re¬ 
gard  to  Journal  bearing  lubrication:  oil  elast¬ 
icity  can  play  as  important  a  role  as  oil  vis¬ 
cosity,  &  shear  rates  much  higher  than  the  value 
10”  S''  (at  present  under  com!  t  tee  consid¬ 
eration)  occur  in  practical  bearing  conditions. 
The  highest  shear  rate  in  the  QM/Shell  studies 
was  5.6  10”s''.  Exxon  studies  (13)  show  that 
even  higher  shear  rates  (above  lO'^s'')  can  occur 
during  acceleration. 

In  the  present  paper,  we  give  information 
about  the  method  used  by  Shell  to  aeasure  Hi  for 
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■ultigrade  oils;  we  describe  recent  atteapts  to 
extend  the  range  of  Ni  measureaent  to  higher 
shear  rates  &  teaperatures  in  order  to  obviate 
extraptolation;  &  we  give  preliminary  results 
which  suggest  that  viscosities  near  2  cP  can  be 
aeasured  up  to  t  =  5  10*8"'.  It  is  better  to 
avoid  extrapolation  if  possible,  particularly 
for  block  copolymer  additives  for  which  morphol¬ 
ogical  changes  occur  at  certain  temperatures. 

ELASTIC  LIQUID  HOLE  PRESSURES 

When  a  liquid  undergoes  shear  flow  between 
parallel  walls,  the  stream  surfaces  are  plane. 
If  one  wall  has  a  small  transverse  slot,  the 
stream  surfaces  near  the  slot  mouth  are  curved 
in  the  sense  illustrated  in  fig. 3.  This  hole¬ 
generated  curvature  coigiined  with  a  shear¬ 
generated  Ni  distribution  gives  rise  to  a 
positive  hole  pressure  P*  :=  Pi -Pa ,  where  Pi 
denotes  the  pressure  exerted  on  the  wall 
opposite  the  slot  center  A  Pa  denotes  the 
pressure  at  the  slot  base.  A  force  balance  on 
the  liquid  contained  within  the  dashed  lines  of 
fig. 3  shows  that  P*  will  increase  with  Ni  . 


FIGURE  3  Hole-induced  curvature  of  the  stream 
surfaces  near  the  mouth  of  a  transverse  slot 
combined  with  a  shear^ flow- induced  Ni  distribu¬ 
tion  generate  a  positive  value  of  P*  =  Pi -Pa . 
Measurements  of  P*  are  here  used  to  determine  Ni 
for  multigrade  oils  at  high  shear  rates. 

In  practice,  there  are  additional  contribu¬ 
tions  to  P*  arising  from  inertial  forces  (be¬ 
cause  the  slot  has  a  finite  width  b)  &  from  the 
response  of  the  flush-mounted  transducer  to  a 
non-uniform  pressure  (because  this  transducer 
has  a  finite  width  (14)  &  may,  in  addition,  be 
imperfectly  aligned  with  the  slot  mouth).  For 
small  values  of  Reynolds  number  Re  A  suffi¬ 
ciently  small  flush  transducer  diaphragm  de¬ 
flections,  it  is  reasonable  to  assume  that  the 
three  contributions  are  additive  A  mutually  in¬ 
dependent,  so  that  we  may  write 

P*  =  P*f  P*ii  +  P*,,  (1) 

where  P*f,  the  flush  transducer  finite  width 
contribution,  may,  for  a  given  transducer  A  die 
geometry,  be  assumed  to  be  determined  by  the 
wall  shear  stress  o  and,  possibly,  the  teaper- 
atuVe;  P**,  the  inertial  contribution,  moy  be 
assumed  to  be  a  product  of  o  A  s(»e  function  of 
Re;  P*» ,  the  "elastic"  contribution,  is  expected 


to  be  a  positive  monotonic  increasing  function 
of  Ni .  Re  is  defined  here  by  the  equation 

Re  =  pbht/(4n),  (2) 

where  p  A  q  denote  the  liquid  density  A  vis¬ 
cosity,  ^  denotes  the  shear  rate  at  the  wall  at 
a  location  unperturbed  by  the  hole,  b  denotes 
the  slot  mouth  dimension  parallel  to  the  main 
flow  direction,  A  h  denotes  the  wall  separation. 

A  "STRESSMETER"  is  a  slit  die  rheometer 
(fig. 4)  with  which  the  hole  pressure  P* ,  the 
shear  stress  o,  A  the  shear  rate  i  can  be  meas¬ 
ured.  Here  we  consider  only  the  case  in  which 
the  hole  cross  section  is  a  narrow  rectangle  (a 
treuisverse  slot);  circular  holes  have  also  been 
used  (11,14,15). 


Ps  Pa 


FIGURE  4  Schematic  diagram  of  a  central  section 
of  a  STRESSMETER  die.  Pressures  Pi,  Pa,  A  Pa  are 
measured  by  means  of  one  flush-mounted  A  two 
hole-mounted  transducers.  The  shear  stress  o  and 
the  first  normal  stress  difference  Ni  are  deter¬ 
mined  from  Pa-Pa  A  Pi -Pa ,  respectively.  The  flow 
rate  Q  is  determined  from  the  rate  of  rise  of 
test  liquid  discharging  into  an  output  cylinder 
(fig. 5).  The  viscosity  is  determined  from  o  A  Q. 

The  "HPBL  theory"  (11,15-17)  relating  P*e  A 
Ni  is  not  as  well  established  as  the  theory  for 
measuring  Ni  in  a  cone-plate  rheometer.  The  HPBL 
theory  for  a  transverse  slot  gives  equations; 

Ni  ^  2nP*e,  (3) 

n  =  d(log  P*»)/d(log  a).  (4) 

Here  A  throughout,  Ni  ,  o,  A  ^  denote  values  at 
the  die  wall  in  a  fully  developed  flow  region 
where  perturbations  from  the  slots’  presence  are 
negligible.  For  the  particular  case  in  which  Ni 
is  proportional  to  o*  over  a  range  which  in¬ 
cludes  zero,  n  =  2  A  Eq.(3)  gives  Ni  =  4P*» ;  if 
b  «  h  A  P*(t  =  0  (a  low  Re  limit),  this  agrees 
with  the  result  of  an  approximate  analysis  (18) 
for  a  particular  form  of  constitutive  equation 
(the  "Second  Order  Fluid");  for  multigrade  oil  A 
many  other  applications,  this  simple  equation  is 
too  restrictive,  A  Eq.(3)  is  more  useful.  For  a 
circular  hole,  Bq.(3)  is  replaced  by  the  similar 
equation  Ni -N2  =  3nP*e  . 

Extensive  tests  of  these  equations,  made  by 
comparing  Ni  A  Ni -Na  values  obtained  from  them 
and  from  measurements  made  with  absolute  rheome¬ 
ters,  have  been  reported  elsewhere  (11,15)  for 
several  polymeric  liquids.  The  curious  result  is 
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that,  within  the  scatter  of  Beasureiient  (usually 
about  ±101!  or  better),  the  range  of  validity  of 
the  HPBL  equations  extends  rather  far  beyond  the 
range  of  validity  of  assumptions  used  in  their 
derivation  (16,17).  The  explanation  of  this 
fortunate  result  is  as  yet  unknown,  but  this 
need  not  prevent  us  from  exploiting  it.  Numer¬ 
ical  simulations  are  also  being  used  to  invest¬ 
igate  the  validity  of  Eq.(3)  for  various 
constitutive  equations  (19,20). 

THB  LODGE  STRESSMBTBR( " >  FOR  HIGH  SHEAR  RATES 

A  schematic  diagram  of  a  central  section  of 
the  STRESSMETER  die  is  given  in  fig. 4.  Located 
far  enough  downstream  from  the  die  entrance,  so 
that  entrance  flow  disturbances  are  negligible, 
are  two  small,  deep,  parallel,  transverse  slots 
of  width  b  <  h,  the  die  height.  Pressures  Pa  & 
Pa  at  the  bases  of  these  slots  are  measured  by 
transducers  Tz  &  Ta .  With  its  effective  center 
located  as  close  as  possible  to  the  Pz  hole  cen¬ 
terline  produced,  a  third  transducer  Ti  is 
mounted  flush  with  the  die  wall  opposite  the 
slot  mouths.  Compressed  nitrogen  or  air  drives 
the  test  liquid  from  a  drive  cylinder  through 
the  die  &  into  an  output  cylinder  (figs. 5, 6)  in 
whose  base  is  a  liquid-filled  hole  connected  to 
a  fourth  transducer  T4  whose  pressure  reading 
time  derivative  gives  the  flow  rate  Q. 

THERMOMETER  PRESSURE  ' 

I - p,.  TRANSDUCERS  T. .  T, .  T,  T.  ! 

.  ‘  OUTPUT  CYUNDER  BLEED  VALXTS  B, ,  Bj ,  B,  I 


SLOT  UNIT 

\ 

\ 


DIE  HEIGHT 
SPACER 


'.'B3  \  Bj 


1  [TEMPERATURE 
I,  DIFFERENCE 
[[THERMOCOUPLE 


FIGURE  5  Schema!  ic  diagram  of  a  STRESSMETER  slot 
unit,  showing  a  flush-mounted  transducer  Ti  , 
slot-mounted  transducers  Tz ,  Ta ,  temperature- 
difference  thermocouple,  downstream  thermometer, 
output  cylinder,  A  flow  rate  transducer  T« .  Typ¬ 
ical  die  dimensions  are  given  in  Table  1. 

Fine-wire  thermocouples  in  contact  with  the 
flowing  test  liquid  are  located  up-  &  downstream 
from  the  die  in  order  to  get  a  rough  measurement 
of  the  teiq>erature  rise  caused  by  the  work  done 
in  forcing  the  liquid  through  the  die;  this  rise 
usually  does  not  exceed  about  1  deg.  C,  which  is 
negligible  for  most  multigrade  oil  applications. 
The  measurement  temperature  is  determined  from 
the  reading  of  a  platinum  resistance  thermometer 
in  the  output  stream  between  the  die  exit  A  the 
output  cylinder  base;  this  is  expected  to  give 
something  close  to  a  cup  mixing  temperature. 
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For  isothermal  shear  flow  through  a  slit 
die  (21),  the  shear  stress  a,  shear  rate  f,  A 
viscosity  n  ore  determined  from  the  equations 

a  =  h(Pz  -  P3)/(2z);  (5) 

^  =  2(2  +  m)Q/(wh2);  (6) 

m  =  d(log  Q)/d(log  a);  (7) 

n  =  ■  0/9.  (8) 

Here,  w  denotes  the  die  width  A  z  denotes  the 
distance  between  slot  centers  measured  parallel 
to  the  main  flow  direction.  Eqs.(6)  A  (7)  embody 
the  well  known  slit-die  analog  of  the 
Weissenberg-Rabinowitsch  correction  for  homo¬ 
geneous  liquids  for  which  n  varies  with  f.  When 
n  is  independent  of  m  =  1  A  Eq.(6)  gives 
f  =  6Q/(wh2). 
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FIGURE  6  Schematic  diagram  of  the  liquid  flow 
system  of  the  high  shear  rate  STRESSMETER. 
Arrows  show  the  flow  path  during  a  run:  with  Z 
the  only  valve  oi>en,  gas  pressure  drives  the 
liquid  fr(»  the  drive  cylinder  through  the 
filter  F  A  slot  unit  into  the  output  cylinder. 
The  other  valves  are  used  for  filling,  bleeding, 
A  emptying.  Heaters  on  the  outer  surfaces  of  an 
oven  A  the  cylinders  are  used  in  conjunction 
with  8  separate  circuits  to  control  temperature. 
The  system  is  surrounded  by  thermal  insulation. 

In  the  ccMKrcial  version  (22)  of  the  high 
shear  rate  STRESSMETER  used  in  the  GM/Shell 
study  (1)  A  in  the  work  reported  here,  specially 
made  diaphragm-capacitance  transducers  are  used 
for  Ti  to  Tw .  T«  is  used  in  a  direct-reading 
mode;  its  signal  is  differentiated  with  respect 
to  time  by  an  analog  circuit.  The  transducer  Ti 
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is  used  as  a  null  indicator  in  a  servo  circuit 
controlling  a  pressiu'e  regulator  which  supplies 
coapressed  air  to  the  electrode  side  of  tne 
transducer  at  a  pressure  which  is  autcwatically 
adjusted  so  as  to  equal  the  pressure  Pi  exerted 
by  the  flowing  liquid.  Sinilar  systeas  are  used 
with  T2  &  Ta.  The  3  liquid  pressures  Pi,  P2 ,  & 
Pa  are  thereby  equilibrated  with  air  pressures. 
Ti ,  T2 ,  &  Ta  do  not  need  to  be  calibrated;  their 
nonlinearity  does  not  aatter;  they  need  stable 
zeros  at  the  working  teaperature.  The  air 
pressures  Pi  &  P2  act  on  the  two  sides  of  a 
high-linear-range  pressure  difference  transducer 
Ts ,  whose  output  gives  Pi  -P2 .  Air  pressures  P2  & 
Pa  are  applied  to  a  second  siailar  transducer 
Ta  ,  whose  output  gives  P2-P3. 

Output  data  read  froa  digital  aeters  are 
typed  into  a  coaputer  for  data  evaluation;  pro- 
graas  have  been  written  which  include  differ¬ 
entiations  required  by  the  "W-R  Equation"  (7)  & 
the  HPBL  Eq. (4) ;  typical  die  diaensions  &  output 
data  are  given  in  Table  1. 

Fig.  6  is  a  scheiMitic  diagraa  of  the  coa- 
plete  liquid  flow  systea.  Filtered  teat  liquid 
is  forced  upwards  into  the  drive  cylinder 
through  the  lower  valve  D,  which  is  then  closed. 
With  valves  R  &  Z  closed,  drive  pressure  applied 
to  the  drive  cylinder  forces  the  liquid  upwards 
through  the  bleed  valves  B2 ,  Ba ,  B4 ,  past  the 
transducers  T2 ,  Ta ,  T4 ,  &  through  the  slots  in 
the  slot  unit.  When  liquid  flowing  into  the 
output  cylinder  is  seen  to  be  bubble-free,  B2-B4 
are  closed  &  Z  is  opened  in  order  to  fill  the 
rest  of  the  die  with  liquid.  During  a  run,  Z  is 
the  only  valve  open. 

After  a  run,  the  drive  pressure  is  reduced 
to  zero  &  the  \)alve  R  is  opened  to  return  the 
liquid  to  the  drive  cylinder.  Density  can  also 
be  aeasured  by  using  the  pressure  transducer  T4 
to  aeasure  the  pressure  generated  by  the  head  of 
liquid  filling  the  tube  connecting  T4  to  the 
base  of  the  output  cylinder.  A  standard  liquid 
of  known  density  is  used  in  a  separate  siailar 
experiaent  at  the  saae  teaperature.  A  coaparison 
of  the  two  pressure  readings  gives  the  required 
density.  An  accuracy  of  the  order  of  IX  or 
better  is  attainable;  this  is  aore  than  adequate 
for  calculating  Re  for  use  in  the  STRESSNETER 
data  analysis. 

The  transducer  T4  can  be  calibrated  in  a 
few  ainutes  in  situ  at  the  working  teaperature. 
The  outer  cylinder  is  filled  with  test  liquid  & 
all  valves  are  closed.  Known  air  pressures 
(aeasured  by  a  Betz  aicroaanoaeter  with  a  range 
of  2.5  kPa)  are  applied  to  the  electrode  side  of 
T4  &  the  P4  readings  recorded.  The  T4  nonlinear¬ 
ity  is  significant,  &  a  linear  second  order 
regression  is  used  to  fit  the  calibration  data; 
about  6  data  pairs  usually  suffice.  For  the  den¬ 
sity  aeasureaent,  T4  can  be  used  solely  as  a 
null  indicator;  the  T4  calibration  is  not  need¬ 
ed.  One-point  calibrations  of  the  pressure  dif¬ 
ference  transducers  Ts,  Te  are  aade  with  a  aer- 
cury  aanoaeter. 


MEASUIffiMENT  OF  P*f  &  P** 

Although  this  section  deals  only  with 
practical  details  of  the  aeasureinent  technique, 
these  are  far  fron  routine  &  are  essential  if 
reliable  aultigrade  oil  data  are  to  be  obtained. 

For  multigrade  oils,  the  average  polymer 
molecular  weight  is  not  very  high  &  the  elastic¬ 
ity  is  rather  small  (compared  with  that  of  other 
polymer  solutions  for  which  Ni  data  are  given  in 
the  literature) ; consequently,  P*f,  P*n  &  P*e 
are  comparable  in  magnitude.  In  order  to  get  P*e 
(and  hence  Ni )  from  the  difference  P*-P*f-P*M, 
it  is  therefore  essential  to  measure  P*f  &  P*n 
accurately.  It  is  reasonable  to  use  the  follow¬ 
ing  equations  for  these  quantities; 

P*f  =  Bior  +  B202  +  ...  +  Bmo«;  (9) 

-P*H/a  =  AiHe  +  AaRe^  +  ...  +  A«Re“.  (10) 

The  coefficients  Bi  nay  vary  with  temper¬ 
ature  and  temperature  history  because  dif¬ 
ferential  thermal  expansion  may  cause  small 
relative  movements  of  electrode  emd  transducer 
case.  If  the  Ti  diaphragm  deflection  gives  a 
negligible  change  in  the  die  height  h,  the  coef¬ 
ficients  Ai  should  be  determined  by  the  detailed 
shape  of  the  T2  slot  edges  alone;  if  these  are 
both  perfectly  sharp,  the  dependence  of  P*ii/a  on 
He  should  agree  with  that  computed  (23)  for  a 
Newtonian  liquid.  In  any  case,  Ai  should  be 
consteuits  for  a  given  slot  unit.  We  assume  that 
P*N  can  be  determined  using  a  Newtonian  liquid 
and  then  used  for  non-Newtonian  liquids  by 
taking  the  value  of  Re  at  the  wall  to  be  equal 
to  that  of  the  Newtonian  liquid. 

Our  present  purpose  is  to  determine  Ai  &  Bi 
experimentally.  If  in  Eq.(lO)  one  changes  the 
independent  variable  from  Re  to  o,  the  coef¬ 
ficients  of  the  powers  of  o  in  the  resulting 
equation  will  depend  on  q-  The  coefficients  Bi 
in  Eq.(9)  do  not  depend  on  n-  In  principle, 
therefore,  one  can  get  Ai  &  Bi  frc*  P*(a)  data 
at  one  temperature  for  two  Newtonian  liquids  of 
different  viscosities.  We  now  consider  what 
values  of  viscosity  q  to  choose. 

ICTBOD  A  -  An  obvious  choice  ("Method  A") 
is  to  take  q  large  enough  so  that  Re  is  suffi¬ 
ciently  small  for  the  inertial  contribution  P*n 
to  be  negligibly  small  over  the  whole  o  range 
required,  so  that  P*f  =  P* .  Having  thus  deter¬ 
mined  P*f,  one  can  then  measure  P*  for  a  Newton¬ 
ian  liquid  of  smaller  q  &  get  I**ii  =  P*-P*f.  We 
have  used  this  method  with  some  success  (11), 
but  we  now  feel  that  it  has  certain  practical 
drawbacks:  the  higher  q  is,  the  more  likely  is 
the  liquid  to  have  a  small  elasticity,  the  more 
difficult  is  the  servo  system  operation,  &  the 
more  time  it  takes  to  change  samples. 

Using  "Bright  Stock"  (q  =  1.3  Pa.s)  at  20C 
to  get  P*f  &  at  higher  temperatures  to  get  P*», 
values  of  At  &  Bi  were  obtained  by  regressions, 
&  were  subsequently  used  to  evaluate  P*»(or)  for 
an  experimental  multigrade  oil  (labelled  "UN20” 
then  t,  ”BC30"  now)  at  5(X:  and  8(K.  A  LODGE 
STRESSMETERC )  (Serial  No.R301)  was  used.  It  was 
found  (11,  fig. 22)  that  the  points  fell  on  a 
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coMKtn  curve  P*»(o)  up  to,  but  not  above. 

Re  =  15.  If  the  tiae-teaperature  superposition 
principle  was  valid  for  this  oil,  these  results 
suggested  that  the  HPBL  eq.  was  valid  up  to,  but 
not  above.  Re  =  15;  this  aeant  that  the  STRESS- 
METER  could  give  values  of  Ni  at  60C  (where 
n  -  22  mPa.s)  for  i  up  to,  but  not  above, 

3  10®s'‘.  The  data  obtained  by  Shell  (1)  with  a 
siailar  STRESSMETER  (Serial  No.  H302)  were 

thought  to  be  subject  to  the  saae  restriction 
Re  <  15;  extrapolation  was  used  to  get  Ni  values 
at  higher  temperatures  &  shear  rates  for  the 
Journal  bearing  studies. 

Our  recent  studies  (again  made  with  Serial 
NO.H301)  cast  doubt  on  the  validity  of  this  re¬ 
striction,  show  no  indication  of  the  existence 
of  a  restriction  up  to  Re  =  99  (the  largest 
value  used  to  date),  &  suggest  that  Ni  can  in 
fact  be  measured  up  to  =  10®s'*  at  lOOC.  Our 
previous  method  of  determining  P*f  and  P*k  in¬ 
volved  the  assumptions  that  P'»  =  0  for  Bright 
Stock  near  20C  &  that  P*f(o)  is  independent  of 
temperature.  Our  new  method  makes  no  such  as¬ 
sumptions. 

METHOD  B  -  Let  P*a  &  P*b  denote  values  of 
P*  measured  using  two  Newtonian  liquids  "a",  "b" 
(preferably  at  one  temperature).  It  is  conven¬ 
ient  to  choose  n  (in  mPa.s)  in  the  range  10-15 
for  one  liquid  &  20-25  for  the  other.  There  is 
a  greater  choice  of  such  liquids  which  ought  to 
be  Newtonian  (i.e.,  P*e  =0),  &  so  additional 
tests  could  be  made  by  comparing  results  ob¬ 
tained  with  low-molecular-weight  liquids  of  var¬ 
ied  chemical  ccxnposition.  If  such  comparisons 
give  consistent  results,  the  assumption  P*e  =  0 
would  be  supported.  Since  the  coefficients  Bi 
should  have  the  same  values  for  the  two  liquids, 
it  follows  from  Eqs.(l),  (2),  (8)-(10)  that,  on 
taking  the  difference  between  values  of  P*(o') 
obtained  for  the  two  liquids,  the  P*f  contribu¬ 
tions  cancel,  leaving  the  result 

■ 

P’b-P**  -  ZAn(bh/4)"{(pa/na2)"-(pt./nb2)'’}o".(ll) 
n  =  1 

It  follows  that  an  m*''  order  linear  regression 
analys  ?  with  X  =  a  &  Y  =  P*b-P*a  will  yield 
the  required  values  of  the  coefficients  An.  One 
can  then  use  Eq.(lO)  for  either  liquid  to 
determine  P*n  &  use  the  P*  data  for  either 
liquid  in  conjunction  with  Eq.(l)  (with  P*e  =  0) 
to  determine  P*f  . 

METHOD  C  -  In  practice,  the  following  pro¬ 
cedure  has  been  found  convenient  with  Method  B. 
We  shall  use  the  term  Method  C  to  mean  measure¬ 
ments  made  according  to  Method  B  &  evaluated  as 
follows.  It  is  convenient  to  use  kPa  as  units. 
Ist  Regression  -  Take  X  =  o,  Y  =  P*a/a,  A  M  =  5 
(say);  determine  Gi .  Gs  ,  where 


{P»b  -  I  GnO"}/or2 

n  =  I 


Z  Dn  Of"  ■  '  . 
n  =  I 


P*./o 


Z  Gn  Of"  ■  *  . 

n  =  1 


2nd  Regression  -  Take  Y  ==  {P*b  -  Z  Gn0f")/of^, 

n  =  I 

X  =  a;  determine  the  coefficients  Di,..D4,  where 


The  required  coefficients  Bi .  Bs  are  then 

given  by  the  equations 

Bi  =  Gi  ,  Bn  =  Gn  +  Dn-l/(l  -  R"-l  ) 

(n  =  2, ...5),  (14) 

where  R  (pbqa^ )/(panb^ ) .  (15) 

P*f  is  now  given  by  Eq.(9)  with  M  =  5. 

One  could  now  determine  Ai  from  the  above 
equations  by  using  the  result 

An  -  Dn/(l-R")/{bhpn/(4n82)}"  (n  =  1,..4),  (16) 

but  it  is  usually  better  to  increase  the  Re 
range  by  combining  P*  data  for  different  liquids 
&  temperatures  &  to  calculate  P*n  from  the  eq. 


P*  -  Z  BnOf". 

n  =  1 


FIGURE  7  Values  of  P*n,  the  inertial  contribu¬ 
tion  to  hole  pressure,  measured  at  three  temper¬ 
atures  using  Method  C,  agree  with  numerical 
simulation  (full  curve)  (23)  up  to  Re  =  30. 

An  illustration  of  results  obtained  by 
Method  C  as  described  up  to  this  point  is  given 
in  fig. 7  for  a  Cannon  standard  viscosity  liquid, 
"S200",  at  50C,  60C,  A  99C.  The  ordinate  is 
P*(i/(aRe),  with  P*m  evaluated  by  means  of 
Eq.(17)  for  the  individual  data  points  at  the 
three  temperatures.  The  full  curve  represents 
results  of  a  numerical  simulation  (23)  with  no 
adjustable  parameter  being  used.  It  is  seen 
that,  although  the  scatter  increases  as  Re  de¬ 
creases,  the  measurements  A  simulation  are  con¬ 
sistent  up  to  Re  -  30;  for  Re  >  30,  the  measure¬ 
ments  lie  above  the  simulation  curve.  (Another 
slot  unit  has  given  data  which  lie  below,  in¬ 
stead  of  above,  the  theoretical  curve.)  Within 
the  scatter,  the  measurements  are  seen  to  be  in¬ 
dependent  of  temperature;  this  was  achieved  by  a 
prior  temperature  cycling  procedure.  Although 
further  experiments  with  intermediate  tempera¬ 
tures  (or  viscosities)  would  be  desirable  for  a 
more  thorough  test  of  temperature  independence, 
the  results  are  evidence  for  the  validity  of  the 
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present  Bethod.  The  difference  between  the  data 
points  &  siaulation  found  for  Re  >  30  aay  be 
due  to  slot  edge  iaperfections  or  to  errors  in  o 
■easureaent  arising  from  aicrobubbles  in  the 
test  liquid,  which  had  not  been  degassed.  Such 
errors  have  been  found,  particularly  at  low 
viscosities  &  high  shear  rates  (see  below). 

For  later  use,  it  is  convenient  to  fit  a 
polynomial  to  experimental  points  such  as  those 
shown  in  fig. 7.  It  is  reeisonable  to  omit  the  few 
outlying  points  &  to  use  the  following  means  to 
obtain  a  reasonable  result  at  low  Re  where  the 
scatter  is  rather  large.  At  low  Re,  the  graph  of 
P*M/(aRe)  V.  log  Re  shows  that,  although  the 
scatter  is  considerable,  the  experimental  points 
are  consistent  with  the  theoretical  curve.  We 
may,  therefore,  use  the  theoretical  low-Re-lim- 
iting  value  -0.035  for  P*M/(orRe)  obtained  by 
several  independent  simulations  (24)  A  force  the 
polynomial  to  give  this  value  as  Re  ->  0.  This 
can  be  achieved  as  follows. 

3rd  Regression  -  Obtain  P*h  values  from  Eq.(17); 
take  X  =  Re°  25,  Y  =  -{P*N/(otHe)  +  0.035}/X;  one 
thus  obtains  a  regrer^sion  equation 

-P*ii/(oRe)  =  0.035  +  azX  + -  +  aeX*, 

X  =  Reo  2s.  (18) 

Use  of  a  power  such  as  0.25  enables  one  to  red¬ 
uce  the  order  of  regression. 
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FIGURE  8  Method  C  values  for  the  flush  trans¬ 
ducer  finite-width  contribution  P*f(o')  are  seen 
to  be  independent  of  temperature. 

Fig.  8  shows  a  final  test  of  this  method 
made  with  the  same  values  of  P*  for  S200  at  50C, 
60C,  &  99C  as  those  used  to  prepare  fig. 7.  The 
ordinate  represents  values  of  P*f/a,  where  P*f 
is  determined  from  P*-F“n,  with  P*h  evaluated  by 
meeuis  of  Bq.(18).  It  is  seen  that  the  data  for 
the  three  temperatures  fall  on  a  common  curve 
(within  the  scatter);  the  full  curve  represents 
the  polynomial  in  Eq.(9),  with  M  =  3.  This  is 
evidence  of  the  validity  of  the  present  method  A 
shows  that,  for  this  transducer  &  die  height, 
P*f(a)  is  independent  of  temperature. 

We  find  that  P*f(o)  can  vary  with  tempera¬ 
ture  after  a  transducer  is  first  assembled  k 
used;  after  a  few  temperature  cycles,  P*f(o) 
usually  settles  dom  to  a  stable,  temperature- 
independent  form  (fig. 8). 


The  final  result  of  using  Method  C  (with  P* 
data  for  two  Newtonian  liquids)  is  the  deter¬ 
mination  of  values  for  the  coefficients  ai  k  Bi 
in  the  regression  equations  (9)  &  (18);  typical 
values  are  given  in  the  footnote  to  Table  1. 

From  P*  data  for  a  non-Newtonian  liquid, 
one  can  then  get  P*e  =  P*-P*f-P*M,  on  using 
Eqs.(9)  k  (18)  for  P*f  k  P*h;  it  is  assumed  that 
the  value  of  n  at  the  die  wall  can  be  used  to 
calculate  Re  for  use  in  Eq.(18).  All  data  given 
below  were  obtained  in  this  way. 

STRBSSMBTER  k  TORSIONAL  BALANCE  RHEOMETER  DATA 

As  a  test  of  our  method  of  measuring  Ni  & 
n,  we  have  compared  STRESSMETER  data  with  data 
obtained  (25)  in  the  UK  on  a  Torsional  Balance 
Rheometer  ("TBR"),  an  absolute  rotational  par¬ 
allel  plate  rhe<»eter  (26,  27)  which  gives  jr ,  a, 
k  N1-N2  for  values  of  f  up  to  2  lO^s"'  .  Samples 
were  drawn  from  the  same  batch  of  a  polymer 
solution  "92"  containing  10.3  wtt  of  polyiso¬ 
butylene  ("Oppanol  B50"  from  BASF;  Mv  =  400,000) 
in  decal in.  Nz  denotes  the  second  normal  stress 
difference  in  shear  flow;  in  order  to  calculate 
Ni  from  Ni -Na  ,  we  used  a  value  0.1  for  -N2/N1  . 
Values  near  0.1  have  been  obtained  (12,  44)  at 
low  shear  stress  for  a  decalin  solution  "Dl" 
with  higher  molecular  weight  polyisobutylene. 
The  "D2"  data  comparison  is  shown  in  fig. 9. 
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FIGURE  9  For  a  polyisobutylene/decal  in  solution 
”D2",  STRESSMETER  k  TBR  data  for  n  &  Ni  agree. 
Over  the  common  t  range  (300  to  20,700  s’ M ,  n 
decreased  from  2  to  0.2  Pa.s,  k  Ni /o  increased 
from  2  to  4. 

It  is  seen  that,  within  the  scatter,  there 
is  agreement  between  the  STRESSMETER  data  k  the 
TBR  data  over  the  limited  coMon  range  of  shear 
stress.  The  small  difference  in  temperatures 
used  makes  no  significant  difference  to  the  com¬ 
parison  because  -d(ln  n)/<)T  is  only  0.01  K’ *  . 

The  agreement  between  the  data  for  ^(o) 
obtained  with  the  two  instruments  suggests  that 
errors  due  to  wall  effects,  flush  transducer 
diaphragm  deflection,  entrance  k  exit  effects,  k 
microbubbles  were  negligible  in  comparison  with 
the  scatter  of  the  data.  The  agreement  between 
the  data  obtained  for  Ni  (o)  is  evidence  of  the 


8 


872043 


validity  of  HPBL  Eqs.  (1),  (3)  &  (4)  but  does 
not,  of  course,  prove  their  validity.  There  may, 
for  example,  be  a  cancellation  of  errors,  &  the 
scatter  in  our  data  may  be  too  great  for  us  to 
detect  errors  such  as  those  arising  from  wall 
effects:  typically,  these  give  differences  of 
about  12\  or  less  in  values  for  viscosity  meas¬ 
ured  in  rotational  and  in  capillary  viscometers 
(28)  when  the  capillary  surfaces  are  specially 
treated. 

STRBSSMETBR  TESTS  FOR  t  >  20,7'H)  s'l 

In  attempting  to  test  our  method  of  meas¬ 
uring  Ni  for  i  >  20,700  s‘ *  (the  TBR  upper 
limit)  we  are  faced  with  the  difficulty  that 
there  is  no  other  rheometer  which  can  give  the 
required  data.  Capillary  jet  thrust  has  been 
investigated  (21,  29,  30),  but  the  difficulties  of 
data  interpretation  have  not  yet  been  over¬ 
come  (21):  Bernd’s  oil  data  (30)  for  f  near 
10®s'*  are  strongly  dependent  on  capillary  di¬ 
ameter  &  are  therefore  unacceptable;  in  addit¬ 
ion,  jet  thrusts  for  Newtonian  liquids  &  multi¬ 
grade  oils  differ  by  only  31,  (31). 

We  must,  therefore,  make  some  assumption 
in  order  to  test  STRESSMETEH  measurements  up  to 
^  -  10®s‘*.  We  shall  use  the  well  known  time- 
temperature  superposition  properties  which  have 
been  established  experimentally  for  a  variety  of 
polymeric  liquids  (32),  notably  those  involving 
flexible  homopolymers  (such  as  polyisobutylene) . 
We  recognize  that  multigrade  oils  containing 
block  copolymers  may  not  give  time-temperature 
superposition.  It  is  reasonable  first  to  test 
the  method  using  polyisobutylene  solutions  and 
others  which  exhibit  time-temperature  super¬ 
position.  If  the  tests  are  positive,  one  can 
then  assume  that  the  STRESSMETEH  method  will 
also  apply  to  block  copolymer  solutions;  one 
should  also  investigate  these  using  TBR  data, 
although  the  t  range  would  be  limited. 

We  choose  the  following  variiifcles: 

dependent:  a,  op ,  or  ,  Opr,  Ni  ,  Nir,  P*»,  P*er; 

independent:  T,  tr .  (19) 

Here,  Op  :=  a  -  (the  polymer  contribution  to 
shear  stress);  n»  denotes  the  solvent  viscosity; 
tr  ^a(T),  where  a(T)  is  a  function  of  the  ab¬ 
solute  temperature  T  alone,  chosen  to  give 
superposition;  8f  Xr  :  =  poToX/(pT)  where  X 
denotes  any  stress  or  pressure  variable  &  po 
denotes  the  value  of  density  p  at  T  =  To . 

Liquids  are  said  to  satisfy  a  time-temper¬ 
ature  superposition  principle  in  ranges 
Ti  <  T  <  Tz ,  tri  <  >r  <  frz  if  a  function  a(T) 
can  be  found  such  that 

dopr/dT  =  0  (20) 

and  dNir/dT  =  0.  (21) 

Such  equations  are  satisfied  by  the  Rotise-Zimm 
theory  for  dilute  polymer  solutions  (33)  when 
certain  model  parameters  (h* ,  N)  are  independent 
of  T  fc  by  Oettinger’s  recent  development  (34)  of 
these  theories  which  gives  9r|/9t  nonzero. 


In  order  to  test  the  validity  of  the  HPBL 
Bq.(3),  we  shall  use  only  the  following  form  of 
time-temperature  superposition  assumption: 

If  ^q.(20)  is  valid,  then  Eq.(^2\)  is  valid.  (22) 


Eq. (20)  &  with  the  equations 
dP*er/9T  =  0 

and  dL*/dT  ^  0, 

9  log  Opr 


where  L* 


9  log  Or 


9  log  tr  /  9  log  >r 

It  follows  from  Eqs. (20),  (22)  -  (25)  that 
9(Ni  -  2nP»»)/9T  =  0, 

9  log  P*»r  /  9  log  Or 


where  n 


9  log 


9  log  tr 


in  agreement  with  Eq.(4). 

Hence,  if  TBR  &  STRESSMETEH  data  show  that 
Ni  =  2nP*e  at  Ti  &  that  the  above  properties 
hold  for  Ti  <  T  <  Tz  St  <  <>  <  .  it  follows 

that  Ni  =  2nP*,  at  Tz ,  and  that  the  HPBL  Eq.(3) 
is  valid  up  to  a  higher  shear  rate  f z  ,  where 


•  1 


(22) 

m 

with 

9 

(23) 

(24) 

(25) 

m  *  1 

(26) 

m  •  1 

(27) 

9  I 


tz  =  tia(Ti)/a(Tz).  (28) 

This  is  the  required  result.  The  only  assuaption 
involved  here  is  Eq.(22),  because  the  range  of 
validity  of  the  remaining  Eqs.  (20),  (23),  S,  (25) 
can  be  found  from  STRESSMETEH  meeisurements. 

If  (as  is  the  case  for  "D2")  the  solvent 
viscosity  n«  is  negligibly  small  in  comparison 
with  the  solution  viscosity  n,  then  opr  -  Or ; 
hence  L*  =  1  St  Eq.  (24)  is  satisfied.  For  multi¬ 
grade  oils,  however,  this  is  not  the  case:  the 
relative  viscosity  nr»i  h/Tb  can  vary  between 
about  1.2  St  2;  for  the  oil  "BC30"  considered  be¬ 
low,  our  measurements  show  that,  while  L*  dif¬ 
fers  significantly  from  1,  it  has  no  significant 
dependence  on  tempterature,  so  that  Eq.(24)  is 
satisfied  within  the  scatter  of  the  data.  This 
result  may  reflect  the  low  ^-dependence  of  n- 

As  an  illustration  of  the  above  procedure 
for  the  case  L*  =  1,  we  use  the  STRESSMETEH  data 
for  "D2''  for  the  range  21C  to  lllC  shown  in 
fig. 10.  The  lower  sets  of  data  points  (repre¬ 
senting  otr(>r)  -  Opr  Or))  dearly  superpose 
(within  the  scatter);  values  of  a(T)  have  been 
chosen  to  give  this  supen>osition.  Thus  Eq.(20) 
is  satisfied.  The  upper  sets  of  data  points 
show,  further,  that  the  same  values  for  a(T) 
also  make  the  data  for  P*rr(i'r)  superpose 
(although  the  scatter  is  somewhat  larger  at  low 
tr),  showing  that  Eq.(23)  is  satisfied.  We  may, 
therefore,  use  the  above  procedure  with 
Ti  =  21C,  Tz  =  me,  li  ^1  -  300  s“ '  ;  we  then 
find,  from  Eq.(28).  that  =  1.1  lO^s  ‘ . 

Using  (22),  we  have  thus  shown  that  the 
STRESSMETEH  gives  valid  Ni  data  for  "02"  up  to 
1.1  10®  S' '  at  11 1C;  here,  n  =  32  mPa.s.  This  is 
based  on  the  comparison  with  TBR  data  at  21C  up 
to  20,700  S'*  (fig. 9);  here,  o  =  3.7  kPa.  The 
STRESSMETEH  data  in  fig. 10  show,  however,  that 
P*»r(tr)  *<  orifr)  are  temperature  independent  up 
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to  t  =  2.9  10®s'i  at  me  (where  n  =  21  aPa.s)  & 
up  to  tr  =  3.9  10®s”*  with  a(T)  =  1  at  lllC. 
This  gives  soae  reason  to  believe  that  the  Ni 
data  are  relied) le  up  to  these  shear  rates. 


LogCPer .  Pa)  +  1 
Log(o'pj',  Pa) 
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FIGURE  10  STRESSNBTER  data  for  log  P*er  (upper 
points)  &  log  Or  (lower  points)  for  the  poly¬ 
isobutylene  solution  "D2"  featured  in  fig. 9  are 
seen  to  satisfy  a  tiae-teaperature  superposition 
principle  over  the  remge  2 1C  to  11 1C.  For  this 
solution,  ot>  =  a  to  within  IX.  The  results 
suggest  that  transverse  slot  hole  pressure  aeas- 
ureaents  coupled  with  the  HPBL  equations  can 
give  values  of  Ni  up  to  f  =  2.9  10®  s'^  (the 
highest  value  for  the  data  shown) . 

To  get  higher  i>  values,  the  polyisobutyl¬ 
ene-decal  in  solution  ”02"  was  diluted  to  give  an 
8.7  wtX  solution  "D2b''.  Siailar  TBR  tests  on 
D2b,  to  be  reported  elsewhere  (35),  support  the 
validity  of  the  STRESSMETBR  aethod  of  deterain- 
ing  Ni  up  to  If  =  1.2  10®s"i  at  99C:  here, 
n  =  10  aPa.s,  Ni/o  =  9.7,  i  Re  =  96;  thus  the 
previous  restriction.  Re  <  15,  is  invalid. 

Ni  DATA  FOR  A  NULTIGRAOB  OIL  "BCSO" 

STRESSMETER  data  for  an  undegassed  aulti- 
grade  oil  (labelled  "BC30"  in  house)  are 
presented  in  figs. 11-14  &  Table  1.  Ni  is  too 
saall  for  TBR  aeasureaent,  so  tiae-teaperature 
superposition  is  again  used  to  test  our  aethod 
of  aeasuring  Ni .  nrfi  varies  froa  1.6  to  1.2,  so 
Hs  is  a  significant  part  of  n  which  varies  froa 
40  aPa.s  at  40C  to  13  aPa.s  at  74C,  the  highest 
teaperature  used.  The  changes  in  n  &  Hp  with  t 
are  about  4X  &  50X. 

Fig.  11  shows  that  P*cr  &  aost  opr  data  tsde- 
en  at  different  teaperatures  superpose.  The  de¬ 
viations  of  the  Opr  data  at  74C  (crosses)  at  the 
highest  fr  are  attributed  to  errors  caused  by 
bubble  trouble  (see  below).  These  outlying 
points  were  oaitted  when  a  regression  was  per- 
foraed  to  express  opr(Fr)  as  a  polynoaial;  a 
regression  was  also  used  to  express  P*fr(5’r)  as 
a  polynoaial.  The  two  polynoaials  were  used  to 
evaluate  derivatives  in  order  to  deteraine  n  by 
aeans  of  Bq.  (27).  The  opr  data  superposition 
shown  by  the  lower  set  of  points  in  fig. 11  shows 
that  Eq. (20)  is  valid,  within  the  scatter. 


We  find  that  o  data  do  not  superpose  when 
the  saae  function  a(T)  is  used.  This  is  iapor- 
temt  here,  because  it  is  not  op  but  a  that  oc¬ 
curs  in  the  HPBL  Eq.(4),  &  so  one  aigbt  think 
that  Ni  data,  evaluated  by  aeans  of  Eq8.(3),  (4) 
would  not  superpose.  However,  the  slopes  of  the 
curves  orifr)  &  orpr(j’r)  have  no  significant  tea- 
perature  dependence  (fig. 12);  in  consequence, 
Bq.(24)  is  satisfied,  n  is  independent  of  teap¬ 
erature,  &  Hi  data  (evaluated  using  Eqs.(3),(4)) 
do  superpose  (fig. 13).  The  validity  of  Bq.(24) 
is  supported  by  the  data  of  fig.  12,  in  which 
log  Or  is  plotted  against  log  opr  with  horizon¬ 
tal  shifts  K(T)  chosen  to  give  superposition. 


FIGURE  11  For  a  aultigrade  oil  at  40C  to  74C, 
STRESSMETER  data  for  P*»r  &  otpr  (the  polyaer 
contribution  to  shear  stress)  can  be  superposed 
by  aeans  of  a  coaaon  function  a(T).  The  opr  de¬ 
viations  at  74C  at  the  higher  fr  are  attributed 
to  errors  caused  by  aicrobubbles  in  the  die. 
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FIGURE  12  Data  for  log  Or  v.  log  opr  can  be 
superposed  by  horizontal  shifts  of  aaount  K  T 
It  follows  that  Eq. (24)  is  valid.  The  scatter  at 
the  high  opr  end  is  attributed  to  errors  due  to 
aicrobubbles  in  the  die. 

The  teaperature  dependence  of  K(T)  tabu 
lated  in  fig. 12  corresponds  to  a  value  of  about 
-0.02  K' *  for  the  relative  viscosity  teaperature 
coefficient  d(ln  nrn )/9T.  This  is  about  20 
tiaes  greater  in  aagnitude  than  values  of  the 


Table  1;  STRESSMETER  data  for  Ni  &  n  for  a  aultigrade  oil  "BC30"  at  40C,  59C,  &  74C 


Run 

Deg. 

DP/SS 

HP/SS 

FHP/SS 

NHP/SS 

EHP/SS 

Re 

Shear 

Shear 

Viscosity 

Ml 

Nl/SS 

n 

[last. 

H-R 

H-R 

No. 

C 

stress 

rate 

iPas 

no,E» 

corr 'n. 

calc 

kPa 

1000/s 

(cP) 

kPa 

10/HPa 

factor 

use'' 

1 

40.0 

220 

0.643 

0.3906 

-.062 

0.320 

1.71 

3.126 

77.70 

40.23 

3.81 

1.22 

1.74 

10.24 

1.021 

yes 

2 

40.0 

221 

0.660 

0.3910 

-.062 

0.332 

1.72 

3.118 

77.67 

40.15 

3.94 

1.26 

1.74 

10.63 

1.021 

yes 

3 

39,9 

225 

0.641 

0.4800 

-.023 

0.189 

0.75 

1.529 

36.04 

42,44 

1.23 

0.81 

1.93 

12.36 

1.026 

yes 

4 

40.1 

219 

0.736 

0.2998 

.097 

0.534 

2.71 

4.711 

119.98 

39.27 

8.93 

1.89 

1.60 

11,33 

1.018 

yes 

5 

40.7 

217 

0.774 

0.2040 

-.129 

0,699 

3.64 

6.344 

lbl.35 

39,32 

14.90 

2.35 

1.49 

11.01 

1.016 

yes 

6 

40.4 

214 

0.740 

0.1025 

-.lb4 

0.802 

4.74 

3.034 

207.13 

38.79 

20,67 

2.57 

1.41 

9.98 

1.014 

yes 

7 

40.6 

213 

0.659 

-.0006 

-.193 

0.353 

5.71 

9.711 

249.92 

38.86 

25.59 

2.64 

1,34 

3.73 

1.013 

yes 

8 

41.4 

223 

0.394 

-.2667 

-.271 

0,931 

8.53 

13.875 

365.36 

37.98 

36.98 

2.67 

1.21 

6.71 

1.011 

yes 

9 

41.3 

212 

0.452 

-.2100 

-.261 

0.924 

3.17 

13.007 

346.03 

37.59 

34.38 

2,68 

1.23 

7.10 

1.011 

yes 

10 

41.0 

212 

0.566 

-.1067 

-.227 

0.900 

6.90 

11.399 

297,80 

38.28 

30.65 

2.69 

1.28 

7.90 

1.012 

yes 

U-R  regression:  In  (Approniate  SR, 

1/lS)  :-l. 

005342E-02  (In  5S)'2  ♦  1 

.085343  (In  SS)  r 

3.104561  ; 

R'2  :  .999909  ; 

SS  In  kPa. 

EHP, Ml  regression: 

In  EHP 

:-. 1592018 

(In  SS)"2 

'  ♦  2.268556  (In  55)  ♦  -2. 

250128  ; 

R'2  :  . 

996571  . 

Density  = 

.8b8  g/il 

1 

59.2 

226 

0.675 

0.4971 

-.082 

0.260 

2.26 

1.213 

56.11 

21.71 

0.64 

0.52 

1.25 

21.31 

1.044 

yes 

2 

59.3 

221 

0.518 

0.4577 

-.142 

0.202 

4.06 

1.931 

94.75 

20.38 

0.87 

0.45 

1.33 

10.48 

1.035 

yes 

3 

59.4 

219 

0,379 

0.3891 

-.233 

0.222 

7,09 

3.152 

159.94 

19.71 

1.73 

0.55 

1.43 

7.05 

1.025 

yes 

4 

59.6 

217 

0.271 

0.2973 

-.329 

0.303 

11.00 

4.754 

244.66 

19.43 

3.85 

0.81 

1.52 

6.38 

1.016 

yes 

5 

60.1 

214 

0.155 

0.1976 

-.402 

0.360 

14.67 

6.451 

329.06 

19.60 

6.54 

1.01 

1.53 

5.53 

1.010 

yes 

6 

bO.5 

212 

0.038 

0.0963 

-.464 

0.406 

13.54 

8.136 

415.48 

19.58 

9.67 

1.19 

1.63 

4.99 

1.006 

yes 

7 

61.0 

207 

-0.136 

-.0178 

-.517 

0.399 

22.73 

9.989 

510.35 

19.57 

12.09 

1.21 

1.67 

4.00 

1.001 

yes 

8 

60.9 

206 

-0.261 

-.1292 

-.553 

0.421 

26.53 

11.752 

597.36 

19.67 

15.40 

1.31 

1.71 

3.59 

0,998 

yes 

9 

61.2 

199 

-0.471 

-.3344 

-.603 

0.466 

34.23 

14.899 

764.58 

19.49 

22,39 

1.50 

1.76 

3.13 

0.993 

yes 

10 

bl.3 

207 

-0.369 

-.2318 

-.585 

0.448 

30.95 

13.342 

687.48 

19.41 

18. 95 

1.42 

1.74 

3,36 

0.996 

yes 

H-R  regression.-  In  (Rpproiiiate  SR, 

1/ls)  :-3. 

018947E-02  (In  SS)*2  »  1 

,143305  (In  SS)  ♦ 

3.768786  ; 

R'2  :  ,9999029  ; 

SS  in  kPa 

EHP, HI  regression: 

In  EHP 

:  ,1210471 

(In  SS)'2 

♦  ,9580033  (In  SS)  ♦  -1. 

,471189  ; 

R‘2  :  . 

991825  . 

Density  -■ 

.855  g/il 

1 

74.2 

22b 

0.270 

0.4903 

-.211 

-0.009 

6.33 

1,341 

99.17 

13.53 

- 

- 

- 

1.044 

yes 

2 

74.2 

229 

0,249 

0.4914 

-.211 

•0.031 

6.32 

1.322 

98.38 

13.44 

- 

- 

- 

- 

1.044 

yes 

3 

74,6 

224 

0.073 

0.4522 

-.319 

-0.061 

10.52 

2  029 

157.25 

12.91 

- 

- 

- 

- 

1.032 

yes 

4 

74.5 

221 

0,034 

0.3910 

-.436 

0.080 

16.67 

3.119 

245.37 

12.71 

0.88 

0.28 

1.73 

2.55 

1.020 

yes 

5 

74,7 

218 

-0.113 

0.2975 

.544 

0  129 

25.55 

4.751 

374.88 

12.67 

2.14 

0.45 

1.69 

2,72 

1.008 

yes 

6 

75.1 

214 

-0.290 

0.1992 

605 

0.116 

34.77 

6.425 

508.57 

12.63 

2.56 

0.40 

1.67 

1.81 

1.000 

yes 

; 

75.4 

209 

-0  356 

0.0909 

-.629 

0,132 

42.32 

8.225 

633  62 

12.88 

5.08 

0.62 

1.65 

2.21 

0.993 

no 

8 

76.1 

207 

-0.462 

-.0169 

- .  636 

0.191 

51.98 

9,974 

774  83 

12.87 

6.39 

0,64 

1.63 

1  91 

0,988 

no 

•) 

76.5 

204 

-0  531 

.1116 

.632 

0,213 

58.22 

11  476 

379.58 

15.05 

8  14 

0,71 

1.62 

1.85 

0.934 

no 

10 

76,3 

199 

-0.644 

-  2667 

.615 

0.237 

68.79 

13.874 

1051.28 

13.20 

10.87 

0.78 

1.60 

1,71 

0.979 

no 

11 

77.0 

195 

■0  b38 

- .  3308 

-  605 

0.248 

73,17 

14,846 

1121  53 

13  24 

12,14 

0.32 

1.60 

1.67 

0.977 

no 

H-R 

regression:  In 

Ihpproiitate  SR. 

1/iS)  :-4, 

183825E 

-02  (In  55)  2  *  1. 

15562’  (In  SS)  ♦ 

4,223731  ; 

R‘2  ,999953  ; 

SS  in  kPa, 

EHr 

,M1  regression: 

In  EHP 

:-4,0877S9E-02  lln  551 '2  ♦ 

1,86902 

(In  55)  • 

-3  437838 

;  R'2 

:  .989662 

Density 

:  .345  g/il. 

:  33.42  iitrofis,  I  -  0523  ci;  t)  licrofl;  v  ;  ,564  ci;  I*  1,'»04  (t'7 
PhP  flush  transducer  contribution  to  hole  pressure  HP,  HHP  inertial  contribution  to  HP,  tHP  ■  HP  -  FHP  -  NHP  D"  ■  drive  pressure. 

HI  1st  norial  stress  difference;  55  :  shear  stress;  [HP  elastic  contribution  to  HP;  [•  :  100  EHP/ISS)  2  10/HPa;  n  ;  d  I09  Hl/d  I09  S5. 

-H«P/SS/Re  ’  0423'4'>f-02  *  .1106126  i  *  1384962  t'2  >  *'3  •  1.994966{'02  X'4  ♦  1 , 318l65( -03  <'5,  uhere  l  :  Re'l,  1  -  .25 

[MP,/S5  5635465  •  5  4010’2[-02  55  ♦  -4  199596E  04  SS'2  »  0  55  3  *  0  55  4;  FHP  (  55  in  kPa  Re  :  rho»b4h45R/(4»etal , 

55  nonlinear  cal  factor  dP/dV5S  1  013593  *  -6.586862E-05  dV55  •  1  3372'>2E-08  dV55'2,  P  in  Pa;  dVSS  in  iV.  norialiied  to  1  at  nO  ci  iercury 

H  R  correction  factor  leiact  shear  rate)/(approiliate  shear  rate);  this  should  be  greater  than  or  equal  to  1 


Hote  1  Due  to  ieasurewnt  errors  in  Runs  13  at  .’4C,  EHP  <  0,  Hi  is  not  calculated  nhen  EHP  i  0 

Hote  2  Runs  ’ll  at  74C  are  not  used  for  the  H-R  correction  because  there  is  a  spurious  viscosity  increase  attributed  to  errors  due  to 

bubble  trouble,  the  oil  uas  not  degassed  prior  to  teasureient 

Hote  3  ipproiiiate  solvent  viscosities  in  iPas:  17.8  at  40C;  11  1  at  59(;  q  35  at  74c 
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siailar  quantity  9(ln  [n])/9'r  obtained  for 
dilute  solutions  of  boaopolyBers  (36,37).  We  do 
not  knoM  the  coaposition  of  "BC30".  It  aay 
contain  a  block  copolymer,  but,  if  so,  then  it 
is  perhaps  surprising  that  time-temperature 
superposition  is  found.  There  is,  however,  some 
uncertainty  about  the  values  of  solvent  vis¬ 
cosity.  For  dilute  solutions,  these  temperature 
coefficients  are  approximately  equal  to  the 
temperature  coefficients  of  the  mean  square  end- 
to-end  distances  of  polymer  molecules  (36). 


which  has  the  dimension  time.  One  interpretation 
of  this  in  molecular  terms  is  given  by  the 
temporary- junction  network  theory  (38),  from 
which,  on  using  the  expressions  given  for  Ni  and 
o,  we  obtain  the  result 

T»  =  <t2>/<t>,  (30) 

where  the  averages  < . . >  are  taken  over  the 
distribution  of  network  junction  lifetimes  t. 
Values  of  x* ,  defined  by  Eq.(29),  are  shown  in 
Figure  14,  and  range  from  8  to  0.4  ps. 


FIGURE  13  STRESSMETER  data  for  Ni  and  Op  for  a 
multigrade  oil  "BC30"  show  time-temperature  sup¬ 
erposition,  except  for  deviations  of  Op  at  74C 
at  high  tr ,  which  are  attributed  to  microbubble 
errors.  The  superposition  is  evidence  of  the  re¬ 
liability  of  these  measurements  up  to  the  high¬ 
est  t  =  1.1  10® s"*;  here,  n  =  13  mPa.s  at  74C. 
The  highest  fr  =  2.7  10®s‘>. 

Using  the  HPBL  Eqs.(3),(4)  to  evaluate  Hi, 
we  find  that,  as  expected  from  figs. 11  &  12,  the 
data  satisfy  the  time-temperature  superposition 
condition  Eq. (21).  Log  Nir(^r)  &  log  apr(tr) 
data  for  40C  to  74C  are  given  in  fig.  13.  The 
superposition  supports  the  view  that  the  STRESS- 
METER  can  give  valid  measurements  of  Ni  for  a 
multigrade  oil  up  to  the  maximum  shear  rate  used 
here;  1.1  10® s‘*  at  74C,  where  q  =  13  mPa.s. 
Using  To  =  74C,  the  reduced  shear  rate  tr 
extends  up  to  2.7  10®s'^.  For  all  these  data, 
Ni/o  <  0.9,  showing  that  the  elasticity  is 
rather  small.  The  greatest  value  of  Re  is  83. 

For  Run  11  at  77C  (Table  1),  it  is  seen 
that  -P*,/P*ii  -  0.248/0.605  =  0.41;  this  is  a 
"figure  of  merit"  for  the  STRESSMETER  method. 
For  a  comparable  quantity  for  the  jet  thrust 
method,  it  is  reasonable  to  use  the  ratio 
(elastic  thruat)/( inertial  thrust);  typical  val¬ 
ues  for  this  ratio  for  multigrade  oils  (31)  do 
not  exceed  about  0.03;  this  is  14  times  smaller 
than  the  above  value  0.41  for  the  STRESSMETER 
method,  which  is  thus  the  better  of  the  two. 

For  analyzing  the  flow  of  oil  in  a  journal 
bearing,  it  is  perhaps  of  interest  to  use  the 
above  data  to  evaluate  the  quantity 

T*  ;=  N>/(2o-t).  (29) 


FIGURE  14  Values  of  log  r*  v.  log  for  the 
multigrade  oil  "BC30"  at  40C  to  74C  evaluated 
from  Stressmeter  data.  r*  :=  Ni/(2orf). 

The  fact  that  the  viscosity  is  not  too 
strongly  dependent  on  shear  rate  suggests  that 
either  the  network  theory  or  the  Rouse-Zin 
theory  might  furnish  a  rough  first  approximation 
to  a  description  of  the  oil;  for  the  Rouse-Zimm 
theory,  in  place  of  Eq.(30)  there  would  be  an 
equation  of  similar  form  involving  sums  over 
discrete  sets  of  relaxation  times  and  their 
squares.  The  choice  between  these  theories  would 
presumably  be  governed  by  the  value  of  c[r|], 
where  c  is  the  concentration  and  [q]  the  intrin¬ 
sic  viscosity,  whose  value  we  do  not  know. 
Oettinger’s  theory  (34)  gives  a  better  approx¬ 
imation,  with  9qp/9t  not  zero. 

VISCOSITY  MEASUREMENT  FOR  y  >  10®s-> 

Commercial  oil  viscometers  have  recently 
been  developed  (39)  which  measure  oil  viscos¬ 
ities  at  150C  and  9  =  10®s'i.  In  the  GM/Shell 
journal  bearing  study  (1),  the  maximum  9  was 
estimated  to  be  5.6  10®s'*.  In  another 
study  (13),  values  above  lO’^s'*  were  found  in  a 
main  bearing  during  acceleration.  It  thus 
appears  that  9  values  higher  than  10®s~‘  should 
be  used  for  current  bearing  A  oil  additive 
research,  whatever  values  may  ultimately  prove 
to  be  best  for  routine  oil  testing. 

He  now  describe  q  measurements  for  9  up  to 
5  10® S'*  for  two  Newtonian  liquids  at  19C: 
degassed  decalin  (2.5  cP)  &  "S6"  (a  9  cP  Cannon 
standard) . 


mi 
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We  use  a  new  "ULTRA-HIGH-SHEAR-RATE  VISCO¬ 
METER";  this  is  a  STRESSMETER  with  a  new  slot 
unit  which  has  a  rigid  wall  in  place  of  the 
flush  transducer;  the  die  length  L,  height  h, 
width  w,  &  slot  width  b  are  reduced  to  about 
2  SB,  54  or  27  pa,  1.6  nn,  &  40  pn,  respec¬ 
tively.  The  reaainder  of  the  STRESSMETER  (fig. 6) 
is  used  unchanged;  it  is  convenient  to  use  an 
annulus  inserted  in  the  output  cylinder  to 
reduce  the  free  liquid  area  in  order  to  Measure 
the  SBaller  Q  values  involved. 


□  CANNON  S6  1st  USE  9JcP 
+  DEGASSED  DECAUN  23cP 


y  ) 


3  4  ' 


FIGURE  15  The  data  points  represent  values  of 
the  viscosity  ratio  rj/oo  at  19C  for  two 
Newtonian  liquids  Measured  in  the  ULTRA-HIGH- 
SHEAR-RATE  VISCOMETER.  ■»  denotes  a  low  shear 
rate  value.  For  decalin,  the  shear  rate  extends 
to  5  lO^s'^.  The  saall  decreases  below  the 
expected  constant  value  line  are  attributed  to 
effects  of  viscous  heating  (fig. 16). 

In  fig. 15,  Measured  values  of  n/Oo  v.  ^  are 
shown  for  the  two  liquids,  oa  denotes  the  low-f 
value.  It  is  seen  that,  above  a  ^  value  which  is 
different  for  the  two  liquids,  n/Oo  decreases 
slightly  as  f  increases.  For  decalin,  the 
greatest  shear  rate  is  5  10^s~^,  t  the  decrease 
there  is  only  8%.  This  result  is  encouraging.  It 
is  expected  that  the  viscosities  of  these 
liquids  will  not  vary  with  shear  rate,  &  it  is 
probable  that  the  observed  SMall  decreases  are 
the  effects  of  non-uniforM  teMperatures  gen¬ 
erated  by  the  work  done  on  the  liquid  in  driving 
it  through  the  die  ("viscous  heating").  To  in¬ 
vestigate  this  possibility,  we  have  replotted 
the  data  as  shotm  in  fig.  16. 

If  (a)  the  observed  decreases  are  due 
entirely  to  effects  of  viscous  heating,  (b)  the 
velocity  &  thersal  fields  are  fully  developed  at 
the  location  of  the  transverse  slots,  (c)  the 
flow  is  laainar,  &  (d)  the  walls  are  isotherMal, 
then  the  values  of  viscosity  ratio  should  agree 
with  those  given  by  the  exact  analysis  of 
Nihoul  (40).  According  to  this  analysis,  the  a 
value  calculated  froH  the  "isotherMal"  Eq.(5) 
should  also  be  correct  for  the  nonisotheraal 
case  &  the  variation  in  n/no  should  be  entirely 
due  to  the  error  involved  in  using  the 


"isotherMal"  Eq.  (6)  (with  b  =  1  for  Newtonian 
liquids)  for  calculating  If  we  now  use  ft  to 
denote  the  isotherMal  value,  given  by  the  eq. 

fi  =  6Q/(wh2)  (isothemal,  Newtonian),  (31) 

then  the  value  t  at  the  wall  in  the  noniso- 
thensal  case  is  given  by  the  equation 

f/fi  =  {1  +  (1  -  2Na)0-5}/2  (Na  <  0.5),  (32) 

where  the  NahMe  nuMber  Na  is  defined  by  the  eq. 

Na  ;=  ph2yo/(36k).  (33) 

Here,  p  :=  -0(ln  n(f,T))/8T  (the  temperature 
coefficient  of  viscosity)  &  k  denotes  the  ther- 
Bal  conductivity.  It  then  follows  that 


n/no  =  2/{l  +  (1  -  2Na)°-5}  (Na  <  0.5).  (34) 


Q  CANNON  S6  1st  USE  9JcP 
+  DEGASSED  DECALIN  25cP 


Log  Na 


FIGURE  16  Measured  values  of  the  viscosity  ra¬ 
tio  n/no  for  two  Newtonian  liquids  plotted 
versus  log  (Nahae  Nuaber)  agree  with  the 
theoretical  curve  for  a  fully  developed  theraal 
field  (isotheraal  walls),  except  for  the  two  S6 
points  at  the  highest  value  of  Na.  The  agreeaent 
suggests  that  the  n/no  decrease  (which  is  only 
8!li  at  Most  for  decalin)  is  entirely  due  to 
viscous  heating  &  that  eui  appropriate  correction 
can  be  aade  by  Means  of  the  theoretical  Eq.(34). 
This  suggests  that  the  ULTRA  HIGH  SHEAR  RATE 
VISCOMETER  can  be  used  to  Measure  viscosities  of 
Newtonian  liquids  in  the  2  to  3  cP  region  up  to 
y  =  5  xl0®8’  * . 

Eq.(34)  is  represented  by  the  full  curve  in 
fig.  16,  in  which  n/no  is  the  ordinate  A  log  Na 
is  the  abscissa.  The  points  represent  single 
values  Measured  for  the  two  liquids.  It  is  seen 
that,  with  the  exception  of  four  points  which 
appear  to  be  the  result  of  randoH  scatter  and 
the  two  S6  points  on  the  right,  the  data  for  the 
two  liquids  fall  on  a  coHMn  curve  which  agrees 
with  the  curve  from  Eq.(34).  The  values  of  the 
constants  for  S6  A  decalin  are  given  in  Table  2. 

The  agreeaent  shown  in  fig. 16  between  the 
Measured  points  for  two  different  liquids  A  dies 
and  the  theoretical  curve  suggests  that  the  de¬ 
crease  in  Measured  q  values  is  indeed  due  to 
effects  of  viscous  heating  up  to  the  greatest  y 
used  for  decalin,  naMely,  5  10*s'>.  Moreover, 
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the  fact  that  two  liquids  having  significantly 
different  values  for  Oi  P,  &  h  (used  with  dies 
of  different  heights)  give  points  on  a  cognon 
theoretical  curve  suggests  that  this  theory  can 
legit iaately  be  used  to  correct  the  aeasured 
data  (interpreted  here  by  ueans  of  isothemal 
equations)  for  the  effects  of  viscous  heating. 
Further  tests  of  this  suggestion  with  other 
liquids  &  choices  of  apparatus  constants  are 
needed. 

The  two  points  for  S6  at  high  Na  which  lie 
above  the  theoretical  curve  in  fig. 16  probably 
reflect  a  departure  from  the  fully  developed 
thermal  field  represented  by  the  full  curve.  A 
thermal  field  is  closer  to  being  fully  developed 
the  smaller  is  the  value  of  the  Graetz  Num¬ 
ber  Gz,  defined  by  the  equation 

Gz  ;=  h3t/(6^),  (35) 

where  a  :=  k/(pCp)  denotes  the  thermal  diffusiv- 
ity.  Thermal  data  are  given  in  Table  2. 


TABLE  2  Thermal  Data,  UHSHV  Teats 


”S6"  Decal in 


Temperature  (C) 

Viscosity  (cP) 

Die  height  (microns) 

(Die  length)/(die  height) 
Maximum  shear  rate  (10*a'i) 
p  (10-2  K-i) 
k  (10-2  n/K/s) 
p  (kg/m2) 

Cp  (103  Nm/kg/K) 
a  (lO--'  m2/a) 

Gz  at  die  exit,  taax 


19 

9.1 

53 

54 
2 

4.1 
13  (?) 

867 
2  (?) 
75 
231 


19 

2.5 

25 

56 

5 

1.7 

11.4 

883 

1.7 

79 

119 


The  larger  Gz  value  for  ”S6"  is  consistent 
with  the  8d)ove  interpretation  of  the  fact  that 
n/lo  lies  above  the  theoretical  curve  for  a 
fully-developed  thermal  field.  It  is,  perhaps, 
surprising  that,  with  Gz  as  large  as  119,  the 
decalin  data  fall  near  this  curve;  the  thermal 
field  must  be  far  from  fully  developed  (41),  but 
it  is  the  velocity  (rather  than  the  thermal) 
field  that  governs  n/n<>  >  the  small  value  of  p 
for  decalin  is  therefore  relevant  here.  Numeri 
cal  simulation  will  be  used  to  investigate  this. 

It  must  be  emphasized  that  these  results, 
encouraging  though  they  are,  have  been  obtained 
for  Newtonian  liquids  only.  For  non  Newtonian 
multigrade  oils,  the  viscosity  varies  signifi 
cantly  (though  not  greatly)  with  shear  rate;  al 
though  this  variation  is  often  no  more  than 
about  30X  up  to  7  =  lO^s  '  ,  the  variation  may 
prove  to  be  greater  at  the  higher  shear  rates 
envisaged  here. 

The  well-known  slit  die  analog  Bqs.(6),  7; 

of  the  Weissenberg-Rabinowitsch  calculation  of 
wall  shear  rate  has  been  derived  on  the  assump 
tion  that  the  test  liquid  is  homogeneous.  This 
Msumption  is  not  valid  when  the  effects  of  vis 
cous  heating  are  significant.  It  will,  there¬ 
fore,  be  necessary  to  see  whether  an  extension 


of  theory  can  be  made  (42)  which  will  allow  one 
to  calculate  wall  shear  rate  from  measured  slit 
die  data  when  the  viscosity  varies  with  shear 
rate  &  the  temperature  is  non-uniform.  The  fact 
that,  for  our  ULTRA-HIGH-SHEAR-RATE  VISCOMETER, 
the  viscous  heating  errors  are  rather  small 
(about  8%  for  decalin  at  5  10®s-*)  gives  reason 
to  hope  that  a  perturbation  analysis,  which 
treats  the  two  corrections  as  independent  of  one 
another,  might  be  sufficiently  accurate  for  ap¬ 
plication  to  multigrade  oils  under  the  condi¬ 
tions  of  journal  bearing  lubrication. 

Other  possible  sources  of  systematic  error 
in  our  viscosity  measurements  include;  o  errors 
arising  from  unequal  values  of  P*n  at  the  two 
slots;  variation  of  n  with  pressure  along  the 
die;  non-parallelism  of  the  die  walls,  further 
tests  are  required  in  order  to  establish  the 
range  of  reliability  of  our  measurements,  but 
our  estimates  show  that  these  errors  should  be 
negligible  in  comparison  with  the  scatter. 

For  the  ULTRA-HIGH-SHEAR-RATE  VISCOMETER, 
the  isothermal  wall  assumption  seems  reasonable 
because  of  the  large  mass  of  metal  which  sur¬ 
rounds  the  very  small  die  region  on  all  sides. 
For  the  STRESSMETER  die,  however,  one  wall  is  a 
thin  stainless  steel  diaphragm  which  will  con¬ 
duct  heat  less  effectively  than  the  solid  stain¬ 
less  steel  wall  opposite  the  flush  transducer. 
The  question  therefore  arises  whether  both  walls 
should  be  assumed  to  be  isothermal  or  whether 
(as  another  extreme  possibility)  one  should  be 
assumed  to  be  isothermal  &  the  other  adiabatic. 
Some  evidence  relevant  to  this  question  is  given 
in  the  next  section. 

MEASUREMENTS  OF  TEMPERATURE  RISE 

Assuming  that  the  die  walls  &  the  liquid 
entering  the  die  have  a  uniform  temperature 
field,  the  effect  of  viscous  heating  is  to  gen¬ 
erate  a  non-uniform  temperature  field  in  the 
flowing  liquid.  At  the  early  stages  of  the  de¬ 
velopment  of  this  field,  the  temperature  profile 
over  a  cross  section  normal  to  the  main  flow 
direction  has  a  maximum  (41)  at  some  point 
between  the  die  center  &  the  die  wall. 

By  moving  the  end  of  the  downstream  thermo¬ 
couple  (fig. 5)  relative  to  the  die  wall,  we  were 
able  to  detect  changes  in  the  temperature  rise 
measured  by  the  temperature  difference  thermo¬ 
couple;  values  recorded  in  this  way  using  Probe 
No. 2  are  shown  in  fig. 17  for  different  values  of 
the  drive  pressure.  Measurement  errors  can 
occur  due  to  conduction  of  heat  along  the 
thermocouple  wires,  heat  generation  by  flow  dis¬ 
turbance  caused  by  the  thermocouple  probe,  A 
the  need  to  locate  the  probe  at  some  distance 
downstream  from  the  die;  the  die  cross  section 
is  too  small  to  allow  one  to  insert  a  probe  in 
the  die,  although  surface  thermocouples  on  the 
die  wall  might  be  used.  The  temperature  profile 
at  the  die  exit  can  be  expected  to  be  convected 
downstream  with  the  liquid  to  a  good  extent. 
Probe  No.l,  made  of  thinner  wire  than  that  used 
for  Probe  no. 2,  does  indeed  give  higher  readings 
(fig. 17)  than  those  given  by  Probe  No.l.  This 
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suggests  that  heat  conduction  along  probe  wires 
was  the  doainant  error.  The  data  also  show  that 
the  aeasureaent  errors  were  large  in  %  terns  & 
that  the  greatest  rise  recorded  was  only  2C. 


FIGURE  17  Measured  temperature  rise  due  to  vis- 
ous  heating  in  a  Newtonian  liquid  in  a  STRESS- 
METER  die.  The  error  due  to  heat  conduction 
along  the  thermocouple  wires  is  illustrated  by 
the  use  of  wires  of  different  thickness  (Probes 
1  &  2).  Two  sets  of  readings  are  shown  with 
Probe  2  in  different  locations  in  a  downstream 
cross  section. 


FIGURE  18  The  curve  represents  the  steady 
increase  in  temperature  rise  due  to  viscous 
heating  as  the  drive  pressure  is  slowly 
increased;  the  probe  position  in  a  cross  section 
downstream  from  the  die  was  previously  adjusted 
to  give  the  greatest  reading.  Tne  filled  circles 
represent  the  isothermal  &  adiabatic  limiting 
values  obtained  by  numerical  simulation  (43). 


Fig. 18  shows  similar  data  for  a  Newtonian 
liquid  (Cannon  "S20")  taken  with  the  downstream 
probe  position  adjusted  so  as  to  give  the 
maximum  reading.  The  continuous  curve  represents 
a  chart  recording  obtained  with  a  slowly-in¬ 
creasing  drive  pressure.  The  two  filled  circles 
represent  values  of  the  theoretical  maximum 
temperature  rise  obtained  (43)  by  numerical 
simulation  for  our  die  dimensions  (L  -  6.8  mi; 
h  -  78  pm)  for  the  two  extreme  cases  in  which 
both  die  walls  are  assumed  to  be  isothermal  or 
adiabatic.  It  is  seen  that  the  measured  curve 
falls  closer  to  the  isothermal  than  to  the 
adiabatic  point.  This  is  consistent  with  our 
expectation  that  the  solid  die  wall  should  be 
approximately  isothermal,  but  the  significant 
measurement  errors  illustrated  in  fig. 17  suggest 
that  limited  weight  should  be  placed  on  the 
quantitative  agreement  between  measured  temp¬ 
erature  rise  &  temperature  rise  calculated  for 
the  isothermal  case.  For  this  reason,  we  believe 
that,  in  order  to  assess  and  correct  for  the  ef¬ 
fects  of  viscous  heating,  more  weight  can  be 
placed  on  the  variation  of  n/qo  with  i  &  its 
comparison  with  the  theoretical  result  (fig. 16). 

BUBBLE  TROUBLE 

The  viscosity  data  shown  in  figs.  15  &  16 
for  the  lower-viscosity  liquid  (decalin)  were 
obtained  after  the  decalin  had  been  thoroughly 
degassed  by  slow  filtration  (a  few  seconds  per 
drop)  into  a  high  vacuum  obtained  with  a  backing 
pump.  Use  of  undegassed  decalin  was  found  to 
give  anomalous  data,  with  the  measured  values  of 
q  passing  through  a  maximum  as  i  was  increased. 
Degassing  by  use  of  a  high  vacuum  without  fil¬ 
tration  was  found  to  be  inadequate,  even  when 
ultrasonics  were  used. 

Fig. 19  shows  data  from  the  first  use  of  an 
undegassed  Cannon  standard  viscosity  liquid  "S3'' 
(3.9  cP  at  20C).  Data  were  taken  with  i  increas¬ 
ing.  The  values  of  q/qo  for  points  at  the  lower 
J  values  (below  1.7  10®s'*)  give  a  standard  de¬ 
viation  equal  to  0.95X  of  the  mean  value.  This 
is  a  measure  of  the  visccmieter  performance.  The 
anomalous  increase  found  at  larger  f  values  is 
not  fully  understood.  We  conjecture  that  somehow 
the  passage  of  test  liquid  through  zones  of 
large  pressure  gradient  causes  agglomeration  of 
microbubbles  (already  present  in  the  liquid)  to 
an  extent  that,  when  the  liquid  reaches  the  die, 
some  of  these  bubbles  have  sizes  coaqiarable  with 
the  die  height.  These  bubbles  then  cause  un¬ 
wanted  changes  in  the  pressure  gradient  near  the 
transverse  slots  &  hence  cause  errors  in  the 
measurement  of  o.  We  have  no  direct  evidence  to 
support  this  conjecture,  except  possibly  the 
observations  (a)  that  thorough  degassing  elimin¬ 
ates,  or  very  substantially  reduces,  the  stream 
of  bubbles  observed  (above  a  certain  flow  rate) 
in  the  liquid  flowing  into  the  output  cylinder  & 
(b)  that,  on  a  succession  of  10  runs  at  one  high 
shear  rate,  the  "S3''  viscosity  anomaly  became 
less  pronounced  (as  one  would  expect  because 
passage  through  the  instrument  removed  some  gas 
from  the  liquid). 
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FIGURE  19  For  f  >  2  10®s'i,  viscosity  ■easure- 
■ents  show  an  anoaalous  increase  attributed  to 
the  presence  of  Bicrobubbles  of  a  size  compar¬ 
able  with  the  die  height.  The  liquid  was  not  de¬ 
gassed  prior  to  use.  Each  point  represents  one 
measurement.  For  the  set  of  7  points  taken  at 
the  lower  values  of  shear  rate  (f  <  1.7  lO^s'*), 
the  standard  deviation  of  n/Oo  is  equal  to  0.95k 
of  the  mean.  rv>  denotes  the  mean  value  of  n  for 
these  points. 

It  is  surprising  that  microbubbles  should 
grow  in  a  region  in  which  the  pressure  is  above 
atmospheric.  Although  negative  pressures  do 
occur  at  the  die  exit  at  the  higher  shear  rates 
(because  of  the  Bernoulli  pressure),  the  pres¬ 
sures  are  estimated  to  be  always  above  atmos¬ 
pheric  at  the  slots,  &  this  has  been  checked  by 
direct  measurement.  Negative  pressures,  with 
concomitant  bubble  formation,  at  the  die  exit 
should  not  cause  errors  in  our  measurements.  It 
appears,  therefore,  that  the  effects  of  pressure 
gradients  in  the  flow  path  leading  to  the  trans¬ 
verse  slots  must  outweigh  the  effect  of  having 
the  pressures  above  atmospheric.  This  is  not  too 
surprising  when  one  considers  that  the  greatest 
drive  pressures  used  (about  500  psi)  are  small 
compared  with  pressures  (estimated  to  be  about 
7000  psi)  required  to  force  all  microbubbles  in¬ 
to  solution;  the  value  7000  psi  is  obtained  from 
the  known  value  of  16000  psi  for  water  (45)  on 
the  assumption  that  the  effect  is  surface-ten¬ 
sion  dominated. 

The  Cannon  liquids  are  standards  when  they 
contain  the  equilibrium  amount  of  dissolved  gas. 
but  not  when  degassed.  For  die  height  measure¬ 
ment,  one  must,  therefore,  use  undegassed  stan¬ 
dard  liquids  only  at  the  lower  f  values  where 
bubble  trouble  is  absent. 

If  our  suggested  explanation  is  correct,  it 
is  to  be  expected  that  rotational  viscosM^ters 
would  not  experience  similar  problems  because 
they  do  not  have  significant  pressure  gradients. 
In  loaded  journal  bearings,  on  the  other  hand, 
there  are  large  pressure  gradients,  &  so  it  is 
natural  to  wonder  whether  bubble  agglomeration 
occurs  with  any  significant  effect  in  the 
positive-pressure  zone.  In  the  negative-pressure 
zone,  generation  of  large  bubbles  (cavitation) 
is  well  known. 
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CONCLUSIONS 

1.  1986  Journal  bearing  oil  film  thickness 
measurements  (1,13)  establish  the  need  for 
obtaining  viscosity  &  shear  (not  bulk)  elastic¬ 
ity  data  for  multigrade  oils  at  shear  rates  up 
to  at  least  6  10®s‘'. 

2.  A  new  conercially  available  slit-die 
rheometer,  the  LODGE  STRESSMETER< " >  FOR  HIGH 
SHEAR  RATES,  is  shown  to  be  capable  of  giving 
such  data  at  shear  rates  up  to  10®8‘‘. 

3.  Data  for  one  oil  give  values  of  an  oil 
time  constant  t*  as  low  as  4  lO  '^s. 

4.  Before  measurement,  oils  should  be  de¬ 
gassed  by  slow  filtration  into  a  high  vacuum  in 
order  to  avoid  errors  due  to  microbubbles. 

5.  First  results  with  a  smaller  die  suggest 
that  viscosities  of  the  order  2  cP  can  be  meas¬ 
ured  at  shear  rates  up  to  5  10®3‘'. 
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NQNBNCUTUIS 

o  Wall  shear  stress  &  shear  rate 

Ni  ,  N2  1st  li  2nd  noraal  stress  differences 
r*  Ni/(2ort) 

Pi.Pz.Pa  Measured  pressures,  fig. 4 
P*  Pi  -  P2  ("hole  pressure") 

P*f  finite  width  contribution  to  P* 

P*»,P*N  Elastic  &  inertial  contributions  to  P* 

n  d{log  P*»>/d(log  o) 

h,  w  Height,  width;  slit  die  cross-section 

z,  L  Slot  separation  &  die  length 

b  Slit  width 

■^1  it  calculated  fro*  iso-themal  equation 

T  Absolute  teaperature 

p,  po  Densities  at  teaperatures  T,  To 

a(T)  Tiae-teaperature  shift  function 

r),  Hs  Solution  &  solvent  viscosities 

nr»i  n/qs  (relative  viscosity) 

no  Value  of  q  at  low  shear  rates 

[q]  Intrinsic  viscosity 

c  Polymer  concentrati' 

Re,  X  pbh^/(4q).  (Re)o  2* 

Op  o  -  qmt 

tr  >8(1)  ("reduced  shear  rate") 

Xr  PoToX/(pT)  for  X  =  a,  Op,Ni,  P*, 

Na,  Gz  ph2^a/(36k),  h3^/(6A) 

k,  a  Thermal  conductivity  &  diffusivity 

Cp  Heat  capacity 

p  -d(log  q(^,T))/dT 

ai ,Ai ,Bi  Coefficients  in  Eqs. ( 18) , ( 10) , (9) 

Di ,Gi  Coefficients  in  Eqs. ( 13) , ( 12) 

P*«,P*b  P*  for  Newtonian  liquids  "a",  "b" 

R  (P6q»2)/(p,qb2) 

Q,  a  Flow  rate,  d(log  Q)/d(log  o) 

L*  Defined  by  Eq.(25) 

MW  Molecular  weight 

Mv  Viscosity-average  MW 

A  :  =  B  A  is  defined  by  the  Eq.  A  =  B 
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